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Abstract 

In this work, we study hole transport in a planar silicon metal-oxide-semiconductor 
based double quantum dot. We demonstrate Pauli spin blockade in the few hole regime 
and map the spin relaxation induced leakage current as a function of inter-dot level 
spacing and magnetic field. With varied inter-dot tunnel coupling we can identify dif¬ 
ferent dominant spin relaxation mechanisms. Applying a strong out-of-plane magnetic 
field causes an avoided singlet-triplet level crossing, from which the heavy hole ^-factor 
^ 0.93, and the strength of spin-orbit interaction ^110 geV, can be obtained. The 
demonstrated strong spin-orbit interaction of heavy hole promises fast local spin ma¬ 
nipulation using only electrical fields, which is of great interest for quantum information 
processing. 
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There is great interest in developing solid-state quantum dot based spin-qubits for quan¬ 
tum information applications.^ Early work on electron spins in III-V semiconductors has 
resulted in the demonstration of single spin isolation,^ manipulation,^ and readout.^ The use 
of semiconductors with strong spin-orbit interactions (SOI) has allowed rapid local control 
of individual spins using only electric fields.^ However, in III-V materials, the ever-present 
hyperfine interaction between electrons and the nuclear spins in the host crystal limits the 
spin dephasing time T^.^ In contrast, electron spin qubits in silicon quantum dots, espe¬ 
cially in isotopically enriched ^^Si, which has almost no nuclear spins, have demonstrated 
extremely long T^.^ However, the well isolated electron spins in silicon, though offering ex¬ 
cellent coherence times, make it difhcnlt to use purely electrical qubit control and also limits 
operation speed. 

Valance-band holes in silicon have even weaker hyperfine interaction, but gain a strong 
SOI due to their p-orbital nature. Heavy hole spins trapped in planar quantum dots for 
quantum information processing have been attracting significant attention since the recent 
theoretical predictions of long spin lifetime,^ reduced hyperfine interaction,^ and all electrical 
spin manipulation.^ Very recently, there have been promising results obtained with single hole 
spins isolated in nanowires. But the properties of holes in nanowires are very different to 
those of holes in devices fabricated from planar p-type MOS (pMOS) structures, silicon- 
on-insulator structures, or heterojunction structures. The most significant difference is 
that in nanowires the hole ground state is the spin-1/2 light hole, whereas in two-dimensional 
planar structures the ground state is the spin-3/2 heavy hole. There have already been both 
theoretical predictions and experimental evidence that hyperfine induced spin relaxation is 
suppressed for heavy holes in 2D systems and quantum dots due to the suppressed hyperfine 
coupling to nuclei,and that spin relaxation can be suppressed by the 2D confinement. 
However, to date there has been very little work on the spin properties of holes in planar 


2 



silicon devices. 


A key technique for implementing spin-to-charge conversion in silicon, and for in¬ 
vestigating spin relaxation mechanisms, is Pauli spin blockade. For two quantum dots 
connected in series, if an excess spin up (down) electron has been trapped in the second dot, 
then ground state tunneling of electrons from the first dot onto the second dot can only occur 
if the additional electron has the antiparallel spin, forming a spin singlet. If the electrons in 
the two dots have parallel spins, forming a spin triplet, current through the double quantum 
dot is blocked by Pauli exclusion. This Pauli spin blockade can be lifted by spin relaxation, 
so measurements of Pauli spin blockade can directly probe spin relaxation mechanisms. 

Here, we present hole transport in a planar Si metal-oxide-semiconductor (MOS) based 
double quantum dot which offers the best of both material systems: Rapid local spin manip¬ 
ulation via the strong SOI of holes; and absence of hyperfine interactions. We demonstrate 
Pauli spin blockade in the few hole regime and identify the dominant spin relaxation mech¬ 
anisms as the inter-dot tunneling coupling is varied. Upon application of a strong magnetic 
field, we observe an anticrossing of the singlet and triplet levels, from which we estimate the 
heavy hole g^-factor and the strength of the SOI. 

Figure 1(a) shows a SEM image of a multilayer gate-electrode Si MOS structure identical 
to the one used in this measurement. The multilayer gate structure enables great flexibility 
of operation from defining a high quality single dot (Supplementary Information, SI) through 
to defining few hole double dots, as shown schematically in Fig. 1(b), studied in this work. 
Lead gates (LG) are strongly negatively biased (Vlg = —4 V) to accumulate holes in source 
and drain reservoirs connecting to the ohmic contacts. Lower gates G1 and G2 define dot 
1 and 2, and control their occupancy. Upper gate G3, in between G1 and G2, is used to 
control the inter-dot tunnel barrier and hence the coupling between dots. The lithographic 
dimensions of dot 1 and 2 are defined by the ^ 30 nm width of gates G1 and G2, and the 
~ 100 width of the adjacent gate LG.^^ 
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Figure 1: (a) Scanning electron microscope image of a device identical to that measured, 
(b) Schematic cross-section of the device when operating in double dot mode. The dot 1 
(2) is defined by lower gate G1 (G2). Upper gate G3 is used to tune the inter-dot tunnel 
barrier. The lead gates (LG) induce source and drain reservoirs. The gate electrodes are 
electrically isolated from each other by thin AlOx layers, (c) Double dot charge stability 
diagrams with Vsd=l mV and Vg3=-1.1 V. The source-drain current (Igd) is measured as a 
function of bias on gates G1 and G2. The black dashed lines joining the triple points are 
superimposed on the diagram, delineating different charge configurations, (d) is a close-up 
of the region marked by red dotted lines in panel (c), showing how the inter-dot tunneling 
can be tuned with gate G3 (indicated in the upper right of each panel). 
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Fig. 1(c) presents charge stability diagrams of the double dot. Although we were unable to 
reach the last hole regime before Isd drops below the background noise, the large variation on 
the separation of the triple points suggests that we are in the few hole limit, where the orbital 
energy level spacing is becoming comparable to the Coulomb charging energy. Fig. 1(d) is 
a zoom-in of the red dashed rectangle in (c), showing how the inter-dot tunnel coupling 
can be tuned with gate G3. At Vgs = —1-1 V (top left of Fig. 1(d)) the tunnel coupling 
is weak, with isolated triple points. Making G3 more negative increases the coupling, so 
that the characteristic honeycomb pattern forms (Vqs = —1.15 V, top right of Fig. 1(d)). 
Finally at Vgs = —1.3 V (bottom right of Fig. 1(d)) the inter-dot tunnel coupling is strong, 
and the states in the two dots hybridize, resulting in almost diagonal lines. The increase in 
inter-dot coupling also causes an increase in the separation between adjacent bias triangles, 
as highlighted by the solid red lines. 

Pauli spin blockade is shown in Figure 2. Figures 2(a) and (b) depict the (m-|-l,n-|-l) 

(m,n+2) inter-dot transition, while Figs. 2(c) and (d) show the (m+l,n+3) (m,n-|-4) 

inter-dot transition, here (m,n) denotes the number of holes in dot 1 and 2. Comparing 
data for positive and negative Vgd we observe current rectification with suppressed transport 
only at positive Vsd [Figs. 2(b,d)]. As expected, we do not observe spin blockade for the 
(m+l,n+2) (m,n+3) inter-dot transition (Supplementary Information S2). Spin blockade 

originates from spin-conserved tunneling; here we label the two relevant highest energy hole 
spins, or more precisely the spin-orbit doublets, by equivalent singlet or triplet states. At 
positive bias, the tunneling of the additional hole onto dot 2 can only happen for tunneling 
between S(l,l) —^ S(0,2) singlet states, and is blocked if the holes are in the triplet state, 
T(l,l). Above a certain inter-dot energy level detuning e, the spin blockade can be lifted 
since T(0,2) can now be accessed from the T(l,l) state. Hence, as shown by the dashed 
trapezoids in Fig. 2(b) and (d), the spin blockade region is dehned by the spacing between 
singlet and triplet states, which, in terms of energy spacing, is the gap between S(0,2) and 
T(0,2). This allows us to extract Ast ~800 /reV for both transitions shown in Figs. 2(b) 
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Figure 2: Bias triangles and Pauli spin blockade. Igd is mapped as a function of bi¬ 
ases on gate G1 and G2 with Vsd=±2.5 mV. Vg 3 is set to -I.IV to give weak inter¬ 
dot coupling, (a)-(b) Hole transport through the double dot involving inter-dot tran¬ 
sition (m+l,n+l) -H- (m,n+2). (c)-(d) Hole transport involving inter-dot transition 

(m+l,n+3) (m,n+4). At positive bias [(b) and (d)], current is suppressed inside the 

singlet-triplet gap (Ast)- The green arrow in (a) depicts the detuning axis direction. 
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and (d). At negative source-drain bias, the spin blockade is absent since the holes are always 
loaded into S(0,2) state and can then exit via S(l,l). 



B (T) B (T) 


Figure 3: Double dot current in the spin blockade regime, (a) and (c) show the double 
dot current as a function of detuning £ and out-of-plane magnetic field B. (a) is the same 
transition as in Fig. 2(b), with weak inter-dot coupling while (c) is the same transition as 
Fig 2(d) with two more holes in dot 2. Vgi(Vg 2 ) is swept to generate the vertical traces 
in (a) [(c)] and then converted to £. The vertical trace are shifted by aligning the £ = 0 
transition to a straight line to eliminate occasional random charge noise offsets, (b) and 
(d) Line cuts at £ = 0 in (a) and (c). The data in (b) and (d) are fitted with models for 
different spin relaxation mechanisms (Supplementary Information, S4 and S5). 


If the detuning is smaller than the singlet-triplet gap £ < Ast, spin blockade can also 
be lifted by spin relaxing mechanisms such as spin-flip co-tunneling,^^ hyperfine interac¬ 
tions [T(l,l) ^ S(l,l)],^^ or SOI [T(l,l) ^ S(0,2)].^^ These will allow a leakage current 
through the double dot when it is in the spin blockade regime. To elucidate the leading 
spin relaxation mechanism, the leakage current is studied as a function of detuning s and 
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out-of-plane magnetic field B. Figure 3(a) maps the leakage current in the spin blockade 
region in Fig. 2(b). Fig. 3(b) shows a line cut around £ = 0 in Fig. 3(a), revealing a peak 
in the source-drain current around B=0. The width of the peak, ^ 400 mT, excludes the 
hyperfine interaction as the dominant spin relaxation mechanism. Instead the peak is well 
described by spin-flip co-tunneling in the regime of weak inter-dot coupling and weak local 
dephasing^^(Supplementary Information S3). From fitting our data to the theory of ref. 26, 
we can extract the hole temperature as ^170 mK and the tunneling rate to the source and 
drain as ^39 //eV. These values are consistent with similar data from almost identical elec¬ 
tron double dot devices, and support the interpretation that the double dot is in the weak 
inter-dot coupling regime. This data also allows us to extract the spin relaxation rate due 
to SOI at higher magnetic field, where the spin-flip co-tunneling is suppressed. The finite 
spin blockade leakage current at higher B {B > 0.5T in Fig. 3(a)) is due to the SOI induced 
relaxation rate F^^, so that we can extract F^^ ^ 2.5 neV, which is consistent with the value 
of ^ 4.6 neV extracted from inter-dot tunneling coupling^^ (Supplementary Information S3). 

To further resolve the influence of SOI on inter-dot spin relaxation, we need stronger 
inter-dot coupling. Increased inter-dot tunneling coupling is achieved by introducing two 
more holes in dot 2, which, by comparing Fig. 2(d) to (b), shows up as further sepa¬ 
rated triple points, higher off-resonance current, and a greater extent of the co-tunneling 
current beyond the bias triangles. The spin blockade leakage current is mapped for the 
(m+l,n+3) ^ (m,n+4) transition in Fig. 3(c). Fig. 3(d) shows a line cut around £ = 0 in 
Fig. 3(c). By increasing the inter-dot tunnel coupling, the peak in Ig^ around 5 = 0 T in 
Fig. 3(b) becomes a dip in Fig. 3(d). The suppression of leakage current at 5 = 0 and the 
rapid lifting of spin blockade with finite field signifies strong SOI induced hybridization of 
the T(l,l) and S(0,2) states. The leakage current dip can be well fitted by a Lorentzian, in 
good agreement with the transport model assuming both strong SOI and strong inter-dot 
coupling^^(Supplementary Information S4). In this regime, the leakage current is limited by 
the spin relaxation rate between (1,1) states, and our fitting yields a relaxation rate of the 
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order 3 MHz, comparable to the values reported for electrons in InAs nanowire double dotA^ 
Since the relaxation rate between (1,1) states is not directly related to the inter-dot tunnel¬ 
ing, the reduced lifetime caused by adding two more holes may be due to an enhancement 
of the SOI, which is highly sensitive to the size of the dot. 




Figure 4: (a) Double dot current as a function of out-of-plane magnetic field and inter-dot 
energy detuning. The map is generated in the same way as Fig. 3(a), but over a larger s 
and B range. Below B =4 T, the bottom two lines are attributed to T_(l,l) ^ S(0,2) and 
T_(l,l) ^ T_(0,2) transitions. To eliminate occasional offsets due to charge noise, traces 
have been shifted vertically by aligning the second resonance transition to a straight line 
(Supplementary Information S6). (b) Energy levels of the double dot in the (0,2) configu¬ 
ration. When the Zeeman energy approaches Ast (Sg^/in^ ^ Ast), the SOI mixes S(0,2) 
with T_(0,2) so that they anticross, (c) Double dot energy diagram as a function of inter¬ 
dot energy detuning £ at finite magentic field. The magnetic field lifts the degeneracy of 
the triplets T+, Tq, and T_. The black dot marks the spin non conserved tunneling with 
the mixing of T_(l,l) and S(0,2) by SOI. The red dot marks the conserved tunneling from 
T_(l,l) ^T_(0,2). 
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The mixing of intra-dot spin states by SOI can be revealed at higher field. Figure 4(a) 
shows a map of the double dot current in the same manner as Fig. 3(a) over a larger range of 
detuning and magnetic field. The lines in the map are resonance current peaks corresponding 
to alignment of energy states between dot 1 and dot 2. Since the inter-dot coupling is weak, 
and the factor is negative,the bottom two lines are ascribed to the T_(l,l) ^ S(0,2) 
and T_(l,l) ^ T_(0,2) transitions as sketched in Figs. 4(b) and (c). The T_(l,l) ^ S(0,2) 
transition, as marked by black dot in Figs. 4(c), moves to higher detuning with increased B 
since T_(l,l) moves to lower energy. The T_(l,l) ^ T_(0,2) transition, as marked by red 
dot in Figs. 4(c), does not move in detuning when B is increased since there is no change in 
spin number. 

At 5 ~ 5 T, we observe an avoided crossing between the bottom two lines in Fig. 4(a), and 
an increase in current amplitude for the lower transition for B> 5 T. As shown in Fig. 4(b), 
T_(0,2) moves closer to S(0,2) with increased B. When the Zeeman energy approaches Ast, 
T_(0,2) anticrosses with S(0,2) due to SOI. If we ascribe the gap in Fig. 4(a) to the SOI 
induced repulsion between S(0,2) and T_(0,2), we can fit the the bottom line to the standard 
anticrossing gap expression. From this fitting, we extract the hole Lande g^-factor as 0.93 
and the single dot spin-orbit gap AgQ of 110 fieY. The extracted heavy hole factor is similar 
to the theoretical prediction of 0.84 for 2D heavy holes,and the values measured in boron 
accepters in Si.^^ Similar values to the extracted AgQ have been reported in InAs^^ and InSb 
nanowires. The spin orbit length can be estimated to be ^ 110 nm,^^ which is considerably 
larger than the size of dot 2. However, we note that some caution should be exercised since 
the two lines identified in Fig. 4(a) do not follow the conventional evolution of a standard 
anticrossing, and there are several resonances at higher £ which cannot be assigned. In 
addition the dependence of these two transitions on in-plane magnetic field (Supplementary 
Information S6) is not well understood, considering the complex spin properties of spin-3/2 
holes. 
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To summarize, we have studied spin dependent transport in a hole double dot system 
based on a highly flexible multi-gate silicon pMOS structure. We observe Pauli spin blockade 
at different charge conflgurations, and observe a leakage current due to spin relaxation. By 
studying the magnetic held dependence of the leakage current, we identify spin-flip cotunnel¬ 
ing as the dominant spin relaxation mechanism at weak inter-dot coupling, while spin-orbit 
driven spin relaxation dominates at stronger inter-dot coupling. Application of a strong 
magnetic field reduces the singlet-triplet splitting, from which we obtain a hole g^-factor con¬ 
sistent with spin-3/2 heavy holes. An avoided crossing is observed, from which we estimate 
the spin-orbit strength. These results give an insight into the properties of heavy hole spins 
in silicon quantum dots and provide a pathway towards spin-based qubits and quantum 
information processing using heavy holes. 

Methods 

Device fabrication. The device studied in this work is fabricated from a high-resistivity 
(p > 10 kfl-cm) natural (001) silicon substrate. The p+ ohmic regions in Fig. 1(b) are 
prepared by solid source boron diffusion at ^ 975 °C with a peak doping density of ^ 
2 X 10^^/cm^. Subsequently, a 5.9 nm gate dielectric (Si02) is grown by dry oxidation in 
the active region of the device. Lower gates G1 and G2, lead gate LG, and upper gate 
G3 are sequentially patterned on the gate dielectric by electron beam lithography, thermal 
evaporation of aluminum, and metal lift-off. Between each patterning, the lower aluminum 
electrodes are oxidized for 10 min at ~ 150 °G in air to form ^ 4 nm AlO^^ for electrical 
isolation. The final stage is forming gas (95% N2/5% H 2 ) annealing to reduce the Si/Si02 
interface disorder and enhance low-temperature performance. 

Experimental setup. Electrical transport measurements are carried out in a dilution 
fridge with a base temperature lower than 30 mK. All signal lines are cold Altered. 
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Supplementary Information: 

Pauli Spin Blockade of Heavy Holes in a Silicon Double Quantum 
Dot 


SI. Single hole transistors with quantum dots defined under different 
gate electrodes 



»Si 18102 ■p+Ohmic Si ■ Si02 ■p+Ohmic 

■ Al 1 AlOx ■ Al 1 AlOx 



Figure SI: (a) and (b) Schematic cross section of the Si MOS multi-gate structure operating as a single 
hole transistor. For (a), a single quantum dot is induced by upper gate G3. Biases on finger gates G1 and 
G2 are set to Vgi=-0.8V and Vg2=-0.83V. For (b), a single quantum dot is defined by G2 with Vgi=-1.4V 
and Vg 3=-2V. Lead gates are set to -4V. (c) and (d) Source-drain bias spectroscopy of the quantum dots 
depicted in panels (a) and (b). For (c), the dot is in many hole regime with weak coupling to source and 
drain reservoirs. For (d), the dot is strongly coupled to the reservoirs with pronounced co-tunneling 
current in Coulomb blockade regions. 

Figure SI shows that with different gate bias configuration, the multi-gate structure can operate 
as a single hole transistor (SHT) with isolated hole quantum dots defined under different gates. 
As shown in Fig. Sl(a), by setting G1 and G2 to depleting mode and G3 to accumulation mode, a 
single quantum dot can be defined under G3 and operate in a similar way as in Ref. SI. Bias 
spectroscopy for this configuration is shown in Fig. Sl(c). The large number of highly periodic 
Coulomb diamonds demonstrates that the dot is in many hole regime and that the Si MOS multi¬ 
gate structure is highly stable. 













Fig. Sl(b) shows that a single quantum dot can be defined under G2 by setting G1 and G3 to 
accumulation mode. As shown in Fig. Sl(d), there is a finite current inside Coulomb diamonds, 
suggesting a high transparency of the source and drain tunneling barriers. This is consistent with 
the geometry of the structure, since the width of the tunneling barrier is defined by the '"4nm 
AlOx gap between Al gate electrodes. Inside the Coulomb diamond in Fig. Sl(d), with Vqs around 
-0.9V, an enhanced conductance can be found when the source-drain bias (Vsd) is higher than 
"'O.GmV. This is corresponding to inelastic co-tunneling via an excited state [S2] and an orbital 
level spacing of "'O.GmeV can be extracted. This value is consistent with the singlet-triplet 
spacing extracted from double dot measurements with a similar plunger gate bias. 


S2. Proof of Pauli spin blockade: double dot charge stability diagram and 
bias triangles 
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Figure S2: (a) Double dot charge stability diagram with upper gate G3 set to -I.IV and source-drain bias 
set to -IV. (b) Detailed gate G1 vs. G2 map showing inter-dot transitions of (m+l,n + l)0(m,n+2), 
(m + l,n+2)0(m,n+3), and (m + l,n+3)0(m,n+4). The source-drain bias is set to -3mV (left) and 3mV 
(right). 


Figure S2(a) shows the double dot charge stability diagram with weak inter-dot coupling. The 
dotted lines joint the triple points. At the top right corner of Fig S2(a) the number of holes is 
small and the source-drain current drops below the background noise so we cannot determine 
the exact hole occupancy. However, even for adjacent charge configurations, there are large 
variations of the size of the parallelograms in this regime, suggesting that we are approaching 
the last few holes in the double dot system. 


To show that the double dot system exhibits Pauli spin blockade we compare the pairs of bias 
triangles with different hole occupation in dot 2. From top to bottom, both maps in Fig. S2(b) 
shows three pairs of bias triangles with consecutive inter-dot transitions of (m+l,n+l)<-^(m,n+2). 









(m+l,n+2)0(m,n+3), and (m+l,n+3)0(m,n+4). Comparing with the negative source-drain bias 
triangles, we see that at positive bias, there is a suppression of the current at the base in the 
bias triangles for the top and bottom pairs in Fig. S2(b). In contrast, by changing the hole 
occupancy on dot 2 by one, there is no suppression of the current for the middle triangle pair. 
Though the spin filling of hole quantum dot can be highly polarized [S3,S4], the alternating 
appearance of suppression current is a characteristic signature of ever-odd spin filling and Pauli 
spin blockade [S5]. 

S3. Lever-arm of the double dot 

Information about the dot size and few hole limit can be obtained from the geometry of the bias 
triangles in Figure 2 in the main text. By the ratio of the plunger gate voltage span to Vsd, the 
gate voltage to dot energy conversion ratio, or lever-arm a, of G1 and G2, can be calculated as 
shown in Table SI. If we follow the semi-classical interpretation [S6], = CGi(G 2 )/Ctotai/ 

which is the ratio of capacitive coupling of G1(G2) to the total environment of the dot 1 (dot 2). 
Since dot 1 and 2 have identical lithographic dimensions, the fact that > a 2 suggests dot 1 
has smaller coupling to the other gates and reservoirs, and hence the dot 1 size is smaller. 
Meanwhile, when two move holes are introduced into dot 2, ^2 exhibits a -^25% decrease, 
implying a either considerable expansion in the size of dot 2, a drastically altered orbital 
wavefunction [S7], or the electrostatic shift of orbital level [S8], all of which further confirms the 
double dot is approaching the last few holes. 

Table SI: Lever-arm a of G1 and G2 for different numbers of holes on the double dot 



(m+l,n+l)^(m,n+2) 

(m+l,n+3)^(m,n+4) 

ai 

0.32 

0.29 

ai 

0.29 

0.21 


S4. Spin flip co-tunneling in spin blockade regime 

The rate equation describing low spin-flip co-tunneling in weak inter-dot coupling regime is [S9]: 
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with the co-tunneling coefficient 
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where e is the electron charge, is the Boltzmann constant, T is the hole temperature, g is the 
Lande g-factor of spin-3/2 heavy holes (0.93 is used, Zeeman splitting is the 

Bohr magneton, 5 is the magnetic field,is the background current offset,/i is the Planck 
constant, (F^^) is the tunneling rate of the source (drain) to dot 1 (2), A is the depth of the 
two-hole level [S9], and U' is the inter-dot charging energy. 



As shown in Figure 2(b), there is no observable change of the current amplitude along the zero 
detuning line, so we assume = F and A U' 1.25meV [Sll]. Then the fitting shown 

in Figure 3(b) yields a hole temperature 7 170 mK and a tunneling rate to source/drain 

F 39 peV. 

The tunneling through the double dot in the absence of spin blockade can be described by the 
equation [S 6 ]: 


lie) = e 



(S3) 


where £ is the double dot energy detuning, F^ (F 3 ) is the rate of a hole entering (exiting) the 
double dot, and t is the inter-dot tunneling rate. Equation S3 can be used to describe the same 
inter-dot transition with negative source-drain bias, i.e. no spin blockade, which gives F^ = F^^ 
and F 3 = F^. By assuming F^ F^^ = F, the fitting result of Eq. (S3) to Isd is plotted in Figure S3. 
The asymmetry of the current peak around 8=0 is due to inelastic tunneling processes at positive 
detuning so only the data points at the rising edge of the peak are used. F 68 peV and 
t 0.56 peV can be extracted from the fitting. 


The tunneling rates to source/drain reservoir F from two different fitting methods agree well 
with each other (fitting to spin flip co-tunneling at positive Vsd gives F 39 peV and fitting to 
elastic inter-dot tunneling at negative Vsd gives F 68 peV). Furthermore, the extracted values 
are consistent with the limits of weak inter-dot tunneling and weak co-tunneling (^f2t < 
k^T, ck^T « F), justifying the use of Eq. (SI) in analyzing the leakage current. 
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Figure S3: (nn,n+2)^(nn+l,n + l) inter-dot transition (non-blockaded, Vsd=-2.5nnV) around zero detuning. 
The black solid line shows fitting result of Eq. (S3) to the £<0 data points. 

While fitting the spin-flip co-tunneling leakage current, a background current= 95 fA can be 
extracted from the fitting. This current continues out to higher magnetic field (with reduced 
amplitude) as shown in Figure 4, and is attributed to spin orbit interaction induce leakage 
current. By a rough estimation assuming that the spin orbit induced singlet-triplet coupling F*^^ 




is the only rate limiting factor in 1^, we can use the relation of to estimate 2.5 

neV. 


According to Ref. [S13], the spin orbit induced relaxation rate can also be estimated based on 
inter-dot coupling as t^/T where F and t are extracted from the fitting of Eq. (S3). This 

alternative analysis yields F*^^ 4.6 neV, confirming the ascription of to spin orbit interaction. 


S5- SOI induced leakage current in the spin blockade regime 

The leakage current at zero detuning due to spin orbit interaction can be described by a 
Lorentzian line shaped dip around zero magnetic field [S13]: 


I = L 


8 ^ , 


(S4) 


with dip width parameter and dip height where F^ei is the spin relaxation rate. 

As show in Fig. 3(d), by fitting the source-drain current in the regime of spin blockade and 
increased inter-dot tunneling coupling to Eq. (S4), we can extract F^ei 3 MHz. 

The background current could be attributed to higher order co-tunneling processes [S6]. 

S6- Double dot current as a function of in-plane magnetic field and inter¬ 
dot energy detuning 

The (0,2) spectrum with in-plane magnetic field is plotted in Figure S4(b). This map differs 
greatly from Figure 4(a) in main text and is not yet fully understood. For in-plane field, the 
resonance lines are considerably wider, which blurs the T.(1,1)^T.(0,2) transition. Extra lines 
within the singlet-triplet gap can be found, which has also been reported in InAs self-assembled 
quantum dot [S14]. However, unlike Ref. S14 where the extra state does not interact with the 
states of interest, a weak anti-crossing can be traced around 2T between the T.(1,1)^S(0,2) and 
the unidentified state. A pronounced crossing around 4T can also be found. However, the 
identification of its origination is challenging without the knowledge of the exact number of 
holes in each dot, and the energy spectrum of holes in 2D systems is not fully understood. 




Figure S4: (a) Double dot current as a function of detuning and in-plane magnetic field. The map was 
taken over a 15-hour period, as gate bias was swept for many different magnetic fields. To eliminate 
the effect of random telegraph noise during the measurement, we fit a parabola to the bottom 
transition in the raw data as marked by the red dashed line in the top of panel (b). The offset applied to 
each sweep relative to the parabola is plotted in the bottom of panel (b), where the two-level 
fluctuator can clearly be seen, (c) Double dot current as a function of detuning and out-of-plane 
magnetic field. The map was the same as Figure 4(a) in the main text, which was taken over a 16-hour 
period. The device stability was much better during this measurement, but charge noise was still 
present as shown in the raw data in the top of panel (d). Since there are several holes in the double dot, 
and the g-factor may differ in the two dots, for simplicity we align the second transition line to a 
horizontal line. The offset applied to each sweep is plotted in the bottom of panel (d). We emphasize 
that this adjustment does not change the energy gap extracted from the data, since the spacing 
between transitions remains unchanged. 
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